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ABSTRACT 

Recent spectroscopic measurements in open clusters younger than the Sun, 
with [Fe/H] > 0, showed that the abundances of neutron-rich elements have 
continued to increase in the Galaxy after the formation of the Sun, roughly 
maintaining a solar-like distribution. Such a trend requires neutron fiuences 
larger than those so far assumed, as these last would have too few neutrons 
per iron seed. We suggest that the observed enhancements can be produced by 
nucleosynthesis in AGB stars of low mass (M < 1.5Mq), if they release neutrons 
from the ^^C(Q;,n)^^0 reaction in reservoirs larger by a factor of 4 than assumed 
in more massive AGBs (M > 1.5Mq). Adopting such a stronger neutron source 
as a contributor to the abundances at the time of formation of the Sun, we show 
that this affects also the solar s-process distribution, so that its main component 
is well reproduced, without the need of assuming ad-hoc primary sources for the 
synthesis of s elements up to A ~ 130, contrary to suggestions from other works. 
The changes in the expected abundances that we find are primarily due to the 
following reasons, i) Enhancing the neutron source increases the efficiency of the 
s process, so that the ensuing stellar yields now mimic the solar distribution at a 
metallicity higher than before ([Fe/H ]> —0.1). ii) The age-metallicity relation 
is rather fiat for several Gyr in that metallicity regime, so that those conditions 
remain stable and the enhanced nuclear yields, which are necessary to maintain 
a solar-like abundance pattern, can dominate the composition of the interstellar 
medium from which subsequent stars are formed. 

Subject headings: Stars: evolution — Stars: AGB and post-AGB — Stars: Carbon — 
Stars: low-mass — Nuclear reactions, nucleosynthesis, abundances 
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Introduction 



Very recently, measurements of s-prqcess elements in you ng Galactic stellar systems 



( iD'Orazi et al 



2009 



Jacobson et al. 



2011 



Maiorca et al. 



2011 



hereafter Paper I) indicated 



that the neutron-rich nuclei Y, Zr, Ba, La and Ce are enhanced as compared to the Sun. 
The enhancement factor is ^ 0.2 dex, and is the same (within uncertainties) for the heavier 
s elements (Ba, La Ce, hereafter hs) and for the lighter ones (Y, Zr, hereafter Is). Hence, a 
roughly solar ratio between the two groups is maintainecjl]. 

The stellar systems involved are open clusters (hereafter OCs) younger than the Sun, 
belonging to the recent evolutionary stages of the Galactic thin disk. In the evolution of 
the disk, over several Gyr, the observed metallicity [Fe/H] (= log (A^(Fe)/A^(H)) — log 
{N (Fe) / N (H.)) q) does not increase appreciably with time (rather, for every epoch the 
observations of [Fe/H] show a considerable scatter). These facts hamper the possibility 
of identifying definite trends i n the abundances of hea vy elements as a function of iron 



( iMashonkina &: Gehren 



2001 



Mashonkina et al 



20071 ). In the cas e of open clusters, 



Magrini et al 



20091), so 



however, one has a reasonably good knowledge of their ages (see e.g. 
that using this information we could show the abundances of neutron-rich elements directly 
as a function of time. This was the key to reveal clearly the mentioned enhancement. 
(Previous studies on fi eld stars had already tentati vely suggested this as a potentially 



effective procedure, see 



Mashonkina fc Gehren 



20011 1. 



Chemical and cheni o-dynamical models of Galactic evolution ( iTravaglio et al. 



2004 



Raiteri et al. 



1999 



19991), based on the currently-accepted scenario for s processing, predict 



a plateau, or even a decrease in the abundances [Xj/Fe] of neutron-rich elements after the 



^The abundance of Sr, in general difficult to obtain, could not be measured because its 
useful lines were outside the wavelength range of the available spectra 
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solar formation; a late increase is certainly not expected. 



The new results can be understood in either of two possible cases. The first possibility 
lays in the late emergence of a purely secondary-like neutron-capture process, i.e. one driven 
by a neutron source not synthesized by the star from pure H, but derived from elements 
inherited at stellar birth. The abundances of s elements thus produced would grow steadily 
in time as compared to iron, becoming appreciable only in advanced epochs. This would be 
the case for the ^^Ne(a,n)^^Mg source, activated in Intermediate Mass Stars (IMS). ^^Ne is 
indeed produced by two a captures on ^'^N, which in its turn derives from p-captures on 
the original CNO nuclei collected from the interstellar medium (ISM). However, the present 
knowledge of s processing in the Sun and in the Galaxy , and especially of the low neutron 



densities at which it must occur f Arlandini et al. 



19991 ). excludes this possibility. Indeed, 



it would imply an isotope mix in the solar composition, for elements sensitive to reaction 
branchings along the s-process path, incom patible with meteoritic abundances (see e.g. 



Kappeler et al 



1990 



Arlandini et al 



19991 ). It would also produ ce an insufficient number 



of neutrons to feed effectively the hs elements (iBusso et al. 



19881 ) a nd would fail to e xplain 



the observed abundance ratios among Sr, Y and Zr in carbon stars (lAbia et al. 



20011 ). 



The only alternative is that the contribution from the ^^C(a,n)^®0 neutron source in 
low mass stars be larger than so far adopted in AGB models. In such a case, however, 
it is important to verify that any model change required to fit the OC data does not 
affect th e many constraint s on A GB nucleosynthesis that have already been accounted for 



(see e.g. 



Busso et al 



1999 



20011 ). In the present paper we aim at demonstrating that, 
for explaining the new observational evidence, it suffices that the ^^C reservoir be larger 
than previously assumed only in very low mass stars (M < 1.5Mq). This would not 



modify previous interpretations of s-enriched AG B stars in the so 



have a higher mass (about 1.5 to 2.5Mq, see e.g. 



Abia et al. 



2002 



ar neighborhoods, which 



Guandalini et al. 



2006 
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Guandalini &: Bussdl2008l ). It must be noted that the standard models used so far for 



them (as well as for the Sun 
to the solar abundances (see 



underproduce 



Travaglio et al 



;he li ght elements Sr, Y and Zr with respect 



2004 . for details). We instead aim at showing 



that our assumptions have the further advantage of avoiding this problem, thus accounting 
completely for the solar main component, i.e. for the abundances of all s elements heavier 
than Kr-Rb. 

This paper is organized as follows. In Section 2 we briefly discuss the changes necessary 
in the models in order to explain the new observational data and we outline our approach. 
In Section 3 we discuss the new stellar yields. In Section 4 we present our main results, 
and substantiate the claim that they also explain the solar s-process distribution. Then, in 
Section 5 we summarize our main conclusions. 



2. Revised models for the s-element enrichment 



In the current literature, one usually finds the same assumptions on the formation 
of the ^^C neutron source, for all LMS between 1 and 3Mq. In particular, the mass of 
the layer where ^^C is formed (the so-called C-pocket) is kept constant in most s-process 
nucleosynthesis models, while the average abundance by mass of ^^C inside the pocket is 
parameterized, to give sufficient neutrons to reproduce the observed constraints. (One has 
actually to recall that there is no ab-initio physical model for the formation of this neutron 
source, but the ob servational data req uire s processing from its activation in AGB stars, as 



discussed, e.g., by 



Gallino et al. 



19981). 



Recently, looking for a better physically-justified scenario, the proton penetration 



into the intershell zone was assumed to be the consequence o 



velocity profile of the convective eddies at the envelope border (jStraniero et al. 



a naturally-decayin 



2006 



2009; 
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Cristallo et al.ll2009[ ). This decay was regulated by a free parameter but 



value (/3 ~ 0.1) was derived for it, independently on the stellar mass (see e.g. 



agam a unique 



Straniero et al. 



20061). 



There is ample space to revise the above assumptions in the attempt to explain the 
s-process abundances in OCs; however, as mentioned, the revisions must preserve the 
interpretation of the abundances in s-process-enriched AG B stars, which have already been 



successfully achieved (IBusso et al. 



2001 



Abia et al. 



20021). As these last are known to have 
masses M > 1.5 M©, we have a large freedom in choosing the parameters for the ^'^C pocket 
in stellar masses lower than this limit. If born today (at solar or higher-than-solar 
metallicity) such low mass stars would barely activate, during the AGB stages, the third 
dredge-up (TDU) necessary to carry He-burning products to the surface. However, the 
efficiency of TDU rapidly increases for decreasing metallicity, so that at progressively lower 
values of [Fe/H] we expect that stars of progressively lower mass can contribute (albeit 
with a limited number of pulses and TDU episodes). One of the scopes of this paper is to 
explore their impact on the Galactic chemical evolution, especially for the recent epochs in 
which OCs were formed. 

We take general pre dictions for mas s loss rates, TDU efficiency and total number of 



pulses from the work by 



Straniero et al. 



(I2OO3I ). This paper provides us with the general 



stellar parameters, which have been often used in works whe re a global picture of AGB 



nu cleosynthesis was need ed, like e.g 
in 



Straniero et al 



Wasserburg et al. 



(120061 ). The stellar models described 



( 120031 ) are also the same used in 



calculations for the s process (ITravaglio et al. 



1999 



previ ous Galactic chemical evolution 



20041 ). By adopting them we can verify 



the effects introduced exclusively by the new proposal here advanced. 

In the present work, we have modeled the formation of the ^^C pocket, for stellar 
masses below l.SM©, by assuming that hydrogen can be injected into the He-rich layers at 
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TDU. The abundance of protons Xh is assumed to exponentially decrease for increasing 
depth (in mass) into the He zone. The ma ximum depth reach ed was set to 4x10~^Mq, i.e. 



4 times larger than currently assumed (cf 



Gallino et al 



1998|). At the restart of shell-H 



burning, the upper part of this region becomes enriched in ^^N, a strong neutron filter 
hampering s-processing. The effective zone where the nitrogen competition is negligible 
turns out to correspond to a total mass of ^'^C available for s-element production M13 = 1.2 
— 2 xlO~^M0 (depending on the details of the proton profile assumed). This means that, 
over the depth of the pocket (4xlO^'^M0), we have an average effective abundance by mass, 
X13, in the range ^ 3 — 5 x 10~^. For mo re massive AGB stars, assumptions very similar to 
those discussed by 



Travaglio et al. 



(120041 ) were instead maintained, with an average mass of 



13 



C burnt per cycle of 3.1xlO"^M0. 

With the above prescriptions we computed the s-process yields (i.e. the contributions 
in solar masses from the AGB stellar winds to the ISM). The set of neutron-capture cross 



sections and of other nuclear inputs was taken from 
further recent improvements; see e.g. iKappeler et al. 



;he KADONIS database and from the 



(120111 ) for a general discussion. 



Our stellar yields were then introduced into an improved version of the Arcetri 



Galactic chemical evo 



Travaglio et al 



(11999 



ution model, which is very similar to the code used and updated by 



20041 ). For the present work we have considerably improved the code. 
This revision included: i) adopting a fine mass resolution (O.OSM©) for the contributing 
AGB stars in the range of interest (1.25< M < 1.5Mq); ii) using a detailed dependence on 
metallicity of the stellar yields (providing them at seven different values of [Fe/H]); and iii) 
estimating the Galactic enrichment using relatively small time steps (25 Myr), in order to 
better describe the final phases of the Galactic chemical enrichment. In particular, sampling 
properly the change of the stellar yields with metallicity is of paramount importance; in 
many previous models the stellar contributions were estimated at few metallicity values 
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(sometimes only two), but such a c ourse grid is not a dequate for following s-element 



abundances. This is clear, e.g., from 



Busso et al. 



(Il999[ ). in particular their Figure 12, where 



it is shown how, for any choice of the ^^C neutron source, the stellar s-element yields do not 
increase monotonically with time and/or metallicity, so that their trends must be sampled 
carefully. Since s-processing is dominated by LMS contributors, also a high resolution in 
time and stellar mass is important, as for them a minimum variation in mass corresponds 



to large variation s in the evolutioii a ry li 



computations by 



Domfnguez et al 



etimes, Tev These last were improved using the 



( 119991 ). which refer to the same FRANEC code from 



which the stellar evolution scenario is taken. 



For the rest, the technique adopted and the main Galactic-model ingredients are very 



close to those of previous works on the sub . 



compared with the ones by 



Travaglio et al. 



ect, so tha. t the results can be immediately 



(1999 



20041 ). In this respect we have adopted. 



for the sake of consistency with preyious works, the solar meteoritic abundances from the 



compilation by 



Anders &: Grevessd (Il989[ ). On the other hand, for internal consistency. 



when referring to differential stellar data, we have used the solar photospheric abundances 
derived in Paper I for Y, Zr, La, and Ce and those from D'Orazi et al. (2009) for Ba. Those 
abundances are indeed on the same scale as those of open clusters, so that no systematic 
offset is introduced. 



3. The stellar yields 

The yields coming from very low mass stars in the model proposed here are different 
than previously assumed, as the larger mass of the ^^C reservoir now adopted affects them 
strongly. The changes concern the total production factors in the winds and their trends 
with metallicity and time. A peculiar characteristic of the cases with a high ^^C mass to 
burn but with a small total stellar mass, hence a limited number of pulses, is that the 
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abundances in the He shell can be sampled by dredge-up before an asymptotic composition 
is reached in He layers. As a consequence, the s-element distribution carried to the surface 
changes greatly from pulse to pulse. In the first 3-4 TDU episodes mainly Is nuclei are 
contributed, while subsequently hs elements become more important. The final abundance 
pattern in the stellar wind is achieved through the averaging operated by convective mixing 
over the continuously changing composition that is dredged-up. In any case, no Population 
I AGB star with M < l.bMc^ is found to experience more than 10-11 TDU episodes (in 



accord with 



Bisterzo et al 



2010h . 



We note however that the predicted yields depend on the efficiency of dredge-up, 



which is an uncertain parameter of any model. In the parametrization by 



Straniero et al. 



(120031 ). eq. (3), it is presented as a function of the metallicity, of the H-exhausted mass 
(Mh) and of its increase during an interpulse period ((5Mh). The values of these last two 
parameters depend crucially on the rates for H-burning reactions, and in particular on that 
for ^^N(p ,7)^^0, the most critica l one. The rece nt changes in th e recom mended value for 
this rate flAdelberger et al.lLoill ). not included in 



Straniero et al. 



(120031 ). certainly affect the 



TDU scenario: in particular, small variations in Mh and SM-^ easily lead to a change in the 
minimum stellar mass contributing at Galactic disk metallicities, moving it inside the range 
1.25 — 1.35M0. Important effects on the amount of dredge-u p are also induced by the set 



Cristallo et al 



20091 ). These 



of opacities used and by the treatment of convection (see e.g. 
are the most typical uncertainties that affect today the AGB inputs to Galactic chemical 
evolution models of s elements. As a consequence, all results must be taken with some 



caution; they can only provide a general guidance and not precise quantitative conclusions. 
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Fig. 1. — Upper panel: The s element contributions from a stellar mass below 1.5Mq, 
expressed as production factors in the stellar winds with respect to the solar composition, 
adopting the extended ^^C pocket described in the text. Lower panel: the same contributions 
as obtained for higher masses with a standard ^^C pocket. 

3.1. Stelleir yields as a function of [Fe/H] 

For values of [Fe/H] relevant for the Galactic disk, our typical yields are shown in 
Figure 1. The upper panel refers to our new assumptions on the ^^C pocket for very low 
masses (M < 1.5Mq), using a ^^C abundance Xis = 4 xlO~^. For comparison, the lower 
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panel shows the yields for stars with M > 1.5Mq, which are obtained with the older 
scenario and the older (smaller) choice for -"^^C mass burnt per cycle. The upper panel, 
in particular, shows the situation for stars between 1.4 and 1.5Mq, which still have a 
rather efficient dredge-up for metallicities up to the solar one. It shows that the stellar 
wind composition, with the new ^^C pocket, is characterized by similar production factors 
for hs and Is elements (as compared to the solar distribution) only at a relatively high 
metallicity ([Fe/H] > —0.2). At lower [Fe/H] values the hs are more effectively produced. 
The distribution is remarkably different than obtained for higher masses (lower panel). 
This is so both for the absolute values and for the [hs/ls] ratios. In particular, the absolute 
yields at high metallicities remain high only for lower mass stars, with the larger ^^C pocket, 
while they are essentially negligible for higher masses, with the smaller, normally assumed 
^^C reservoir. 

The relevance of the trends shown in Figure 1 can be understood with reference to 
the age- metallicity relation of our model, illustrated in Figure 2. It is derived a priori, 
from the reproduction of the main species, like oxygen, the a-rich elements and iron and 
its behavior is typical in all modern Galactic evolution models. For our purposes. Figure 
2 shows two important things: i) the SFR has remained rather high throughout the last 
10 Gyr, after its peak in earlier epochs; ii) during the last 7.5 Gyr of the Galactic history 
the metallicity [Fe/H] has grown very little: over this whole, long era, it has remained 
within ±0.2 dex from solar. When considered together. Figures 1 and 2 tell us that most 
of the stars in the Galactic disk were born with high metallicity; moreover, those having 
mass below 1.5Mq (dominating the Initial Mass Function) generated, in their AGB phases, 
s- process yields close to those at the extreme right of Figure 1 (upper panel). The larger 
yields characterizing lower metallicity stars in Figure 1, with element ratios very far from 
the solar distribution, played only a marginal role in setting the composition of recent 
Galactic stellar systems, because too few stars were born in those conditions. These few 
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stars could synthesize effectively only elements whose overabundances in the stellar winds 
were extremely high. It is known that this is the case with lead, for which low- metallicity 



AGB stars provide the so-called strong s-process component (iGallino et al. 
about 50% of 208pb. 



19981 ) ■ producing 




10 5 

Galactic Age 

Fig. 2. — The continuous curve shows the age-metallicity relation obtained by our model, 
i.e. the history of the iron increase in time produced by the Galactic evolution calculation. 
The dotted line reports the Star Formation Rate (SFR), i.e. the mass going into stars over 
the whole Galactic disk at each epoch, in units of Mq/jt. The values of this function are 
divided by a factor of ten, in order to fit it in the ordinate scale. 
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3.2. The Stelleir yields as a function of time 
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H 1 h 
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-6 -4 -2 

Time since formation of the Sun (Gyr) 



Fig. 3. — The same enhancement factors in the steUar winds shown in Figure 1, plotted as a 
function time. The zero of the scale indicates the formation of the Sun. Upper panel: during 
more than 3 Gyr, before the solar formation, the production factors of s elements in the winds 
of very low-mass AGB stars have remained very similar. Lower panel: more massive stars, 
with a smaller ^^C pocket, could not play a significant role at high metaUicities, because of 
their negligible yields. 



In order to make the above statements clearer we can actually transform Figure 1 into 
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a plot as a function of time, using the age-metallicity relation of Figure 2 and setting the 
zero of the scale at the moment of the solar formation. This plot is presented in Figure 3. 
As one can see, the small region at the right side of Figure 1 (for [Fe/H] values larger than 
about —0.2), where the yields of Is and hs nuclei are about similar (wit hin 0.2 dex), is now 



shown to start typically 3 Gyr before the solar formation. According to 



Dommguez et al 



(119991 ). such an interval allows for the complete evolution of stars in the range 1.4 — l.SM©, 
whose lifetime, in this metallicity interval, ranges from 2.5 to 3 Gyr. As mentioned, at the 
relevant [Fe/H] values these stars did have TDU; hence they could affect the composition 
of the Sun and of stars born subsequently with about equal production factors for Is and 
hs elements, thus clarifying how these two groups of nuclei can increase their abundances 
in OCs, while maintaining their ratios relatively constant. We emphasize that, as Figures 1 
and 3 show, at the times and metallicities mentioned, more massive stars had essentially no 
role on s-processing, because their yields were far too low. 

The above figures make clear why the role of LMS in the chemical evolution of 
s elements can change remarkably relative to previous estimates, if a more extended 
^^C pocket is assumed. 



3.3. The yields of a typical single AGB model 

The Is and hs contributions to the Galactic enrichment derive obviously from many 
stars of different mass and metallicity. Nevertheless, we can identify a typical s-element 
nucleosynthesis progenitor, providing yields for Is and hs elements whose mutual ratios are 
very close to those of the subsequent Galactic average (and of the solar distribution). With 
our large ^'^C pocket, this turns out to be an AGB star of 1.4— 1.45M0, with [Fe/H] ~ —0.1. 
By normalizing to the abundances of s-only nuclei produced by such a model star, we can 
infer s-process percentages in its winds for the elements analyzed here. Using this estimate 
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for explaining the solar abundances we obtain 86%, 84%, 80%, 85%, 70% and 76%, for Sr. 
Y, Zr, Ba, La and Ce, respectively (see also later. Table 1). From this, assuming that the 
complement to one of the percentage abundances of the heavier elements (Ba, La and Ce) 
derives from the r process, we infer r-process components of 15%, 30% an d 24% for Ba, 
La an d Ce, respectively. These estimates are slightly higher than those by 



Bisterzo et al 



torn . 



We now present the results on the abundance evolution of open clusters induced by our 
new hypotheses and then briefly discuss the consequences for the solar distribution. 



4. The resulting s-process composition of the Galactic disk 

4.1. Interpreting the observations of open clusters 

In order to compute model expectations for the heavy elements in OCs we must 
necessarily include in the calculations estimates for their production in other astrophysical 
environments, outside the range of LMS to which our models refer. As already mentioned 
in Subsection 13. 3[ for Ba, La and Ce one can infer a rough estimate for their r-process 
contributions from the yields of our model star best reproducing the solar main component. 

For lighter nuclei, the yields from IMS were estimated assuming that no ^^C burning 
occurs in them: the neutrons derive only from the activation of the ^^Ne(a,n)^^Mg reaction 
in thermal pulses. They induce a secondary- like neutron-capture process, hence the yields 
vary as a function of the metallicity. The total IMS production of s elements is very 
uncertain, as these stars evolve almost entirely in the infrared during the thermally-pulsing 
phase, so that observational constraints are minimal and no trustable calibration of their 
mass loss rates, suitable for their use in stellar models, exists. As a consequence of this, the 
number of TDU episodes they undergo is unknown. While, from nucleosynthesis models in 
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thermal pulses, we can derive the s-element distributions in their winds pulse after pulse, 
not knowing the number of these last that can occur makes their global contribution to the 
chemical evolution of the Galaxy a rather free parameter, which must be calibrated ad-hoc. 
The ^^Ne(a,n)^^Mg source is also activated during core-He and shell- C burning phases of 



massive stars (providing the so-called weak component, see e.g. 



Busso &: Gallino 



1985h. 



'or es timating the corresponding yield s we modified previous estimates by iRaiteri et al. 



( 19931 ) following the recent analysis by 



Pignatari et al. 



(|2010| ). performed with updated 



cross sections for a solar-metallicity 25M0 star. The production in stellar winds was then 
scaled at different metallicities assuming a secondary-like behavior. 

Despite the fact that the uncertainties affecting s processing in MS and IMS are very 
large, they are not critical for our purposes, as their joint contributions to Sr, Y and Zr 
do not account for more than about 10% o f the total solar inventory for these species 



(ISerminato et al 



2009 



Bisterzo et al. 



2010l ). Concerning the r-process (or, in general, the 



primary) components of Sr, Y and Zr, in the presence of huge uncertainties we estimated 

leir low values, t hese r ather arbitrary choices are not 



Serminato et al. 



( I2OO9I ) gave for them values between 5 



them to be around 5%. Due to t 
critical in any respect. Note that 
and 11%. 

In order to compute the Galactic chemical evolution of the five observed elements, 
we proceeded in the following way: i) we derived metallicity-dependent stellar yields for 
the weak component in massive stars and for the IMS (3-8Mq). ii) We then assumed that 
the r-process components of the studied demeri ts are p roduced at the lowest mass end 
of type II Supernovae, as in iRaiteri et al.l (119991 ) and in 



Travaglio et al. 



mm . iii) With 



the above ingredients, we followed the chemical evolution of the Galaxy with our model, 
fixing the uncertain value of the total mass contribution from IMS so that, at the epoch 
of the Sun's formation and at the corresponding Galacto-centric radius, the sum of all 
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non-LMS contributions gives the complement to one of our "best average" AGB model (see 
Subsection I3.3p . This means that the percentages form astrophysical sources different from 
IMS must be of 16% for Y, 20% for Zr, 15% for Ba, 30% for La and 24% for Ce (see also 
Table 1, column 4, Case C). iv) After this calibration, we introduced the yields of LMS 
AGB stars, thus computing the total abundance of the studied species as a function of time 
in the Galactic disk, up to the recent eras when OCs were formed. 

With the above assumptions, the Galactic chemical evolution model produces, for the 
five measured elements, the evolutionary trends shown in Figure 4, where the OC data are 
shown as empty circles; each of them represents the ave rage of the abundan ce over all the 



Maiorca et al. 



201l|). The data for 



stars analyzed in the cluster (usually five; for details, see 
field stars are indicated with crosses of the same size as the errorbar and are taken from the 
SAGA database. In the left panels, adopting the common way of looking at the chemical 
evolution of heavy elements, we plotted the observations as a function of the iron abundance 
or metallicity, [Fe/H], while in the right panels we showed them as a function of time. 
As already mentioned in Section 1, the plots as a function of metallicity give a confused 
pattern. This is so because the large scatter in the iron abundances at any given age in the 
last few Gyr of Galactic life does not make [Fe/H] a good evolut ionary indicator fo r these 



stages. This finding is in agreement with previous indications by 



Mashonkina et al. 



f!2007h . 



Only when the evolution in time is properly reconstructed (right panels) the recent 
growth of s-elements is unveiled. This is however easily done only for the clusters. Their 
age is kno wn from a varie t y of m ethods: our specific choices, discussed in Paper I, were 



taken from 



Magrini et al. 



( 12009[ ) . For clusters olde r than 0.5 Gyr, th ey made use of the 



morphological age indicator 6V (jPhelps et al. 



19941 ). as calibrated by 



Salaris et al. 



(l2004h : 



for younger clusters they used the most recent age estimates available, like those obtained 
by the lithium depletion boundary method. For field stars, similar criteria do not exist and 



- 18 - 



a precise value of their age is usually not available, so that observations referring to them 
(present in the left panel) cannot be plotted as a function of time. 

The s-element abundances in OCs displayed in Figure 4 were measured differentially 
with respect to the solar photosphere, using the same method of analysis and the same line 
list (see. Paper I for details). For this reason, the data were normalized to the photospheric 
solar abundances. These are, in some cases, slightly different from the meteoritic ones 
(that were used in our "average model" of the solar main component and in deriving the 
s-process and r-process contributions to the chemical evolution). This explains why, for 
some elements, the model curves pass above the solar position: it is only an appearance, 
deriving from the normalization. 

In each panel of Figure 4, the shaded areas represent predictions from our calculations 
and aim at giving an idea of the model sensitivity to the most crucial parameters, i.e. the 
abundance of ^^C burnt and the lower mass limit for TDU operation. The shaded ranges 
take into account both a variation by ±25% around the average value of our ^^C abundance, 
X-12, = (4 ± l)xlO~^, and the uncertainty on the lowest mass contributing to the process, as 
discussed in Section 3. In the right panels, two more evolutionary curves are plotted: the 
long dash-dotted line shows the evolutionary trend obtained by excluding the masses below 
I.5M0; the short dash-dotted one represents instead a case in which they are included, 
but using for them the same extension of the ^^C pocket adopted for higher mass stars. It 
should be clear, from the right panels of the figure, that the extra production we assumed 
from low mass AGB stars offers a means to account for the recent observations of young 
OCs. 

It is relevant to mention that the spectroscopic observations in clusters, which are 
our original constraints, do not offer us indications with an accuracy better than about 
±0.1 dex; within these limits we can say that the chemical evolution model, despite 
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all its uncertainties, suffices to account for the data, namely for a global increase in the 
abundances that nevertheless maintains a solar-like distribution. This is the main point we 
want to stress: s-element abundances in Population I stars can be explained, within their 
uncertainties, by suitably fixing the parameters controlling the neutron source in low mass 
stars, without the need of further ingredients. 

4.2. LMS yields and the soleir distribution of heavy elements 

As a further test, we ran a few cases in which we included in our Galactic model only 
the yields of Figures 1 and 3, excluding any further, uncertain contribution from other 
sources (IMS, MS and the r process). This allowed us to derive an estimate of the role 
exerted solely by LMS in determining the composition in heavy elements of the recent 
Galaxy and of the Sun. This role is illustrated in Table 1. There, we put in evidence the 
uncertainties in the models mentioned before, by showing two possibilities (called cases A 
and B, presented in columns 2 and 3), where we include (or, respectively, exclude) in the 
evolutionary calculations the lowest mass bin from 1.25 to 1.3M0 (see the discussion in 
Section 3 for this). 

As a comparison, we also show as model C, presented in column 4, the s-process 
percentages obtained in the yields of our already-mentioned single AGB star, of 
1.4— 1.45M0 and with [Fe/H] ~ —0.1, best approximating the Galactic average and the 
solar main component. 

Despite the necessary cautions related to the models, the exercise shown in Table 1 
(cases A and B) makes clear that pure s processing in LMS with a suitable choice of the 

parameters for the ^"^C pocket can account for typically 80% (sometimes more) of all the 
solar s-element abundances. This is so for both Is and hs nuclei, without distinction. As 
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Table 1 

Percentage Contributions to Solar Heavy Elements from LMS 



Element 


This work (A) 


This work (B) 


This work (C) 


Literature 


Strontium 


90 


89 


86 


55(1) 


Yttrium 


81 


80 


84 


62(1) 


Zirconium 


80 


78 


80 


55(2) 


Barium 


89 


86 


85 


84(2) 


Lanthanum 


72 


70 


70 


70(2) 


Cerium 


79 


78 


76 


81(2) 



Cases (A), (B): chemical yields from LMS only, with/without the lowest mass bin 
Case (C): co mposition of the wind s o f a single AGB star, l A5Mp,, [Fc/H] = —0.1 



Travaglio et al 



(|2004l ): 



,- 2 



Bisterzo et al 



POIOI ). single star 



Table 2 

Percentage contribution to solar s-only nuclei 



Isotope 


Case A 


Case B 


Arlandini et al.(1999) 


iooRu 


95 


93 


95 


"OCd 


97 


95 


97 


i24xe 


94 


92 


91 



a comparis on, in Table 1 (column 5) we report the percentages of Is elements attributed 



to LMS by 



Travaglio et al 



(|2004| ). without the ^^C pocket enhancement. These authors 
underlined that 20-30% of the abundances were not accounted for, even after including the 
contributions from IMS, massive stars and the r-process. In our new scenario, depicted in 
Table 1 (columns 2 and 3) the LMS production of the elements considered is increased 
just by the amount that was missing (20 - 30%), so that the solar abundances of Sr, Y 
and Zr are now satisfactorily reproduced. We cannot therefore confirm the need, suggested 
by other authors, for an unknown primary process, sometimes c alled LEPP, ot Light er 
Element Primary Process, contributing to these elements (see e.g 
details on this proposal.) 



Travaglio et al. 



2004, for 



In order to substantiate our statements about s processing, Table 2 reports our LMS 
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contributions to a few s-only nuclei of the solar distribution with A < 130. The results 
are presented for the cases A and B defined above. Again, considering the 10% errorbar in 
solar abundances, plus all the model uncertainties, it seems that the solar concentrations of 
s-only nuclei can be rather satisfactorily explained. In particular. Table 2 shows that our 
stellar yields, integrated through the chemical evolution model up to the formation of the 
Sun, can account for essentially the sam e percentages that cou ld be previously obtained 



only from average single-star models (see 



Arlandini et al 



19991 ) 



Note that we consider Tables 1 and 2 only as indications that the solar s-process 
elements can indeed be explained by pure s-processing. It would be completely outside our 
scopes to search for a really detailed fit to the solar abundances through a full chemical 
evolution procedure. 



5. Discussion and Conclusions 

In the previous sections we presented a tentative interpretation of the increase in 
heavy element abundances after the Sun's formation, as revealed by recent spectroscopic 
observations of young open clusters. The increase affects elements that are predominantly 
produced by slow neutron captures and we showed that their trend in recent stellar systems 
can be understood as a result of the contributions from very low mass, long-living stars 
(M < 1.5Mq, Tgt, > 2 Gyr), if they activate the main neutron source more effectively than 
previously suspected. 

A consequence of the above suggestion is that also for the Sun the s-process main 
component can be naturally accounted for, without the need of any extra contribution. An 
additional conclusion that might be obtained from our procedure concerns the behavior 
of predominantly-s elements at low metallicity, in comparison to those that are mainly 
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produced by the r-process, like e.g. Eu. In our cases A and B of Table 1 we found, for 
the s-process component of Eu, the values of 7.5 and 6.5%, respectively. From this, and 
making the guess that the complement to one of the abundances for Ba and Eu is due 
to the r-process, one might deduce a ratio [Ba/Eu],. ranging from —0.8 to —0.9. This is 
roughly compatible with the average values measured in r-process-enri ched metal-poor 



stars, excluding th e largely scattered measurements at [Fe/H] < —2.5 (,Suda et al. 



2008 



Sneden et al 



20081 ). However, due to the mentioned uncertainties, we believe that such an 
extrapolation to the r-process is not really robust (although it is usually accepted in the 
literature). The stellar and Galactic parameters cannot simply be fixed accurately enough 
to yield a real confidence in such estimates. 

Understanding how our scenario is produced is now easy, in the light of the results 
shown in Figures 1 to 3. What is crucial in this context is the evolution in time: most 
parent AGB stars were born in the last 7.5 Gyr, when the metallicity is high and suitable, 
with the larger ^^C pocket we suggest, to produce Is and hs elements in similar proportions 
(similar to those of case C in Table 1). 

We notice that, regardless of the specific model we propose here, the same observational 
data offer a picture of the heavy element evolution in the Galactic disk that is incompatible 
with the hypothesis that fast primary processes strongly contribute to Is elements. Indeed, 
the fact that the abundances of Y, Zr, Ba, La and Ce grow with respect to iron (and do 
so very late in the Galactic history) excludes a primary process, the products of which 
would follow the trend of iron itself. Moreover, the suggestion that a supplementary process 
modifies selectively the lighter neutron-capture nuclei is contradicted by the fact that in 
young stellar systems Y and Zr show the same trend as Ba, La and Ce. Furthermore, below 
A = 130 there are a few s-only isotopes, shielded from fast decays, whose production in 
adequate concentrations requires necessarily the action of slow neutron captures, not of fast 
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processes. 



It is known that the LEPP has been invoked also (and mainly) for extremely metal 
poor stars. Here the classical idea that fast neutron c aptures (thos e of the r process) 



account for the whole abundances of heavy el ements ( 



Truran 



19811) was subsequently 



20081 ). First suggestions on 



integrated and modified by new observations (iSneden et al. 
the existence of different stellar sites for producing r-process elements of different atomic 
mass number were a dvanced from an analys is of the record of live radioactive nuclei in the 



Early Solar Nebula fjWasserburg et al. 



19961 ). A large observational basis on metal-poor 



and extremely metal-poor stars now supports the fact that, while above Te-Xe the r-process 



does appear to be rather universal, the abundances of 



dispe rse and more difficult to interpret (lAoki et al. 



2005 



ower mass e 



Suda et al. 



ements are more 



2008 



Roederer et al. 



20101 ). It was noticed t hat certain heavy nuclei a ppear to be produced in progenitors 



forming very little iron ( IQian fc Wasserburg 



20071 ): it was also suggested that elements at 



the neutron magic number = 50 might derive from a combination of charged-particle 
reactions an d neutron-captures in high-entropy neutrino winds a s sociated to core-colla pse 



supernovae (IQian fc Wasserburg 



Thielemann et al. 



2007 



Wasserburg fc Qian 



2009 



Roederer et al. 



2010; 



20101 ). It is in any case beyond doubts that lighter and heavier r-process 



nuclei should come from different sources or from differerit sup erpositions of astrophysical 



conditions (see e.g. 



Farouqi et al. 



2010 



Thielemann et al 



20111 ) 



The above confused situation is today the object of a lively debate. Unfortunately, 
the data and models discussed in this paper cannot offer any clue for the understanding of 
extremely metal-poor stars. We notice that, also in that case, attempts at explaining the 
complex observational scenario with suitable combinations of existing mechanisms, without 



invoking new ad hoc processes, already exist (|Qian fc Wasserburg 



20091: 



Farouqi et al. 



2010 



Banerjee et al 



20111 ) 



2008 



Wasserburg fc Qian 
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We conclude by mentioning that the extended ^^C pocket assumed in this paper would 
require the existence of very efficient mixing episodes in evolved stars of m ass below l.SMfT, . 



A sirn ilar conclusion was already inferred from very different constraints (jPalmerini et al. 



2011a 



b|); in particular, it seems to be necessary for explaining the isotopic admixture of 



oxygen in presolar oxide grains. The transport mechanisms most commo nly-adopted so 
far (ro tationally-driven shear, or thermohaline mixing) would not suffice (iPalmerini et al. 



2011al ). From this point of view, the apparent need of forming extended ^^C reservoirs in 
the short time interval of a TDU episode, might be another indication that other mixing 
processes must be at pl ay. Phenomena o f this kind were recent ly suggested to be driven 



by magnetic buoyancy (IBusso et al. 



2007 



Nordhaus et al. 



20081 ): this hypothesis and other 



mechanisms must now be examined in order to see if it can offer a general interpretation of 
non-convective mixing in evolved stars. 
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constructive criticisms, which greatly helped in shaping the paper as it appears now. Many 
colleagues provided us with their advise; this is in particular true for C. Abia, A. Chieffi, S. 
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Fig. 4. — The Galactic enrichment in Is and hs elements from our chemical evolution model, as compared 
to observed data. The solar photospheric abundances are shown as heavy dots. The left pane ls show the 
abundances as a function of [Fe/H]; the data for single stars are from the SAGA database by 



Suda et al 



(j2008[ ): their errorbar is not always available, b ut certainly non- negligible. As an example, the estimates 



deducted from Table 3 in 



Simmerer et al 



(|2004[) for La yield an uncertainty of 0.1 dex. Even assuming this 
value (rather optimistically) as typical for all the measured elements, as done in the plot, we see that the 
emerging situation is very confused and a definite trend for all the species studied does not emerge. Simply, 
[Fe/H] is not an a dequate parameter to describe the evolution of s elements. The data for DCs are taken from 



Paper I and from 



D'Orazi et al 



([20091). Each open symbol represents a different cluster; several among them 
have metallicities higher than solar. In the right panels the abundances are instead plotted as a function 
of time: here, only the clusters can be displayed (see text). While the data are limited in number, they 
appear to show clear regularity, without a wide scatter. The shaded area, indicating our model predictions, 
reproduces the solar and clusters' measurements within the uncertainties of the observations and of the 
models (see text for explanations). In the right panels, dash-dotted lines show the trends obtained without 
the extension of the ^'^C pocket suggested in this paper, either including (short dashes) or excluding (long 
dashes) masses below 1.5Mq. 
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